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ABSTRACT 

We report the results of abundance analyses of new samples of stars with planets 
and stars without detected planets. We employ these data to compare abundance- 
condensation temperature trends in both samples. We find that stars with planets 
have more negative trends. In addition, the more metal-rich stars with planets display 
the most negative trends. These results confirm and extend the findings of Ramirez 
et al. (2009) and Melendez et al. (2009), who restricted their studies to solar analogs. 
We also show that the differences between the solar photospheric and CI meteoritic 
abundances correlate with condensation temperature. 



1 INTRODUCTION 



Melendez et al. ( 2009 1 were the first to detect a significant 



correlation between elemental abundances and condensation 
temperature (Tc) in a sample Sun-like stars; in a follow-up 
study, Ramirez et al. ( 2009) confirmed their findings. In par- 
ticular, they found solar twins/analogs to be enhanced in re- 
fractory elements relative to the Sun. Although this was the 
first time such a trend had been found, it had been searched 
for unsuccessfully multiple times before within the context of 
the self-enrichment hypothesis I Huang et al.|2005 Gonzalez 
|2006[ |Ecuvillon et al.|2006[ ) . |Ramirez et al.| ( |2009| speculate 
that the trends are due to planet formation processes. 



In order to test the findings of Ramirez et al. ( 2009 1 



and also to expand on their analysis over a broader range in 
Tcff , we revisit the topic of abundance trends with Tc among 
stars with planets (SWPs). We do so with a new method of 
analysis we introduced in Gonzalez ( 2008 1 and new stellar 



samples, which we described in Gonzalez et al. (20101 



Our paper is organized as follows. In Section 2 we de- 
scribe the new SWP and comparison star samples. We em- 
ploy them in Section 3 to search for evidence of abundance 
trends with Tc; we also examine the recent abundance data 
of [Neves et al.|p009[ ) and the Solar System abundances. We 
discuss our results in Section 4 and present our conclusions 
in Section 5. 



2 DESCRIPTION OF DATA 



In Gonzalez et al. ( 2010 1 we described our most recent SWP 



and comparison stars samples and presented the results of 
our new fine abundance analyses. We compared the Li abun- 
dances of the SWPs to comparison stars using the method 
of comparison introduced in |Gonzalez| ( |200^ and confirmed 
that SWPs near the solar temperature have lower Li abun- 



dances than similar stars not known to harbor Doppler de- 
tectable planets. We deferred reporting the results of our 
analyses of other elements in these samples until the present 
paper. For detailed descriptions of the observations, data re- 
duction and abundance analyses, the reader is referred to 



Gonzalez et al. (20101 



Our initial SWP and comparison stars samples contain 
85 and 59 stars, respectively. We list their abundances in 
Table 1 (provided online). In order to have the best chance 
of detecting the subtle trends between abundance (as [X/H]) 
and Tc, we only calculated the abundances of those elements 
represented by 2 or more quality absorption lines in our 
spectra: Na (2), Al (2), Si (6), Ca (2), Sc (2), Ti (5), V (6), 
Cr (3), Fe (53), Co (3), Ni (6). We employ these data in our 
analyses below. 



3 SEARCHING FOR TRENDS WITH Tc 
3.1 New Data 



In Gonzalez ( 2008 1 we introduced a new index, Ai, which is 
a measure of the distance between two stars in Tcft-[Fe/H]- 
log g-Mv space. In that study we calculated a weighted av- 
erage Li abundance difference between a given SWP and all 
the stars in the comparison sample, with (Ai)~^ serving as 
the weight. We also applied this method (with the addition 
of corrections for bias) in Gonzalez et al. (20101. We employ 



essentially the same analysis method here, only that we are 
comparing [X/HJ-Tc slopes rather than Li abundances. 

The Tc values are from |Lodders et al.| ( [2009| . The ele- 
ments we measured span the Tc range 958 K (Na) to 1659 
(Ca, Sc). This is about the same range of Tc that Ramirez 
|et al.| ( |2009[ ) reported finding an [X/H]-Tc trend. 

We further restricted the samples in the present anal- 
ysis by eliminating stars that do not satisfy the following 
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Figure 1. Weighted-average [X/H]-Tc slope differences between 
SWPs and comparison stars using the samples in the present 
study. The open circle with the dot represents the solar value. 
The error bars for the typical SWP slope difference are shown in 
the lower left. 



criteria: 5400 < Tefi < 6300 K, uncertainty in Tefi < 70 
K, available Hipparcos parallax. After this culling, the SWP 
and comparison stars samples contain 65 and 56 stars, re- 
spectively. 

For each star in our two samples, we calculated the slope 
between abundance (as [X/H]) and Tc using standard lin- 
ear least-squares; the slope values are listed in the online ta- 
ble. The average uncertainty in the slope is about ±0.00007 
dex and is the same for the SWPs and the compar- 
ison stars. We plot the weighted average [X/H]-Tc slope 
differences between SWPs and comparison stars in Figure 
1. Overall, SWPs appear to display more negative [X/H]- 
Tc slopes than the comparison stars, the primary exception 
being some SWPs near the solar temperature. 

When interpreting the datum corresponding to the Sun 
in each of Figure 1 and Figure 2 (and additional figures in 
the following sections), it is important to keep the following 
points in mind. First, although the [X/H]-Tc slope for the 
Sun is identically zero by definition, the weighted difference 
slope does not fall on 0.0 in the figures. The fact that it is 
negative tells us that the Sun's slope is more negative than 
the slopes of similar comparison stars. Second, the uncer- 
tainty in the [X/H]-Tc slope for the Sun is also zero, but the 
difference slope plotted in the figures will have some uncer- 
tainty associated with it due to uncertainties in the slopes of 
the comparison stars. The effective uncertainty of the Sun's 
difference slope should be considerably less than ±0.00007 
dex K~^. This implies its difference slope is significantly dif- 
ferent from zero. 

In [Gonzalez et al.| ( |2010[ ) we introduced an additional 
step in the analysis method of Gonzalez] (j5008 ) to correct a 
systematic bias inherent in the method. We have applied the 
same correction procedure to the present data and present 
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Figure 2. Same as Figure 1 but corrected for bias. See text for 
details. 



(20101 for a detailed description of the correction method). 
The basic trends evident in Figure 1 are not changed. 

The average SWP [X/H]-Tc slope value from the data 
in Figure 2 is (-5.6x10"^ ± 11 (s.d.) 1.4 (s.e.m.)) x lO"'^ dex 
K~^. A simple count yields 46 SWPs with negative slopes 
and 19 with positive slopes. If we assume that slope dif- 
ference values between -0.00007 and ±0.00007 dex K"^ are 
not significantly different from zero and count only the val- 
ues below and above this range, then there are 33 SWPs 
with negative slopes and only 8 with positive slopes. 

Several stars in our samples can be considered solar 
analogs, which we define as having Tcff, log g, [Fe/H], and 
Mv values within ± 100 K, 0.2 dex, 0.1 dex, and 0.2 mag- 
nitude, respectively, of the solar values. Two SWPs (HIP 
15527, 80902) and 9 comparison stars satisfy our criteria. 
We plot the average relative abundances, as [X/Fe], against 
Tc for the SWP and comparison solar analogs in Figure 3. 

The slopes of least-squares fits to the trends in Figure 
3 are (4.0 ± 3.3) x 10"^ and (9.6 ± 6.1) x 10"^ dex K"^ for 
the comparison and SWP solar analogs, respectively. This 
shows that SWPs have an average slope 2.4 times that of the 
comparison stars. This is very close to the ratio of 2.2 found 
by |Melendez et al.| |2009^. While our slope values cannot be 
directly compared to theirs (they included a larger range in 
Tc), it should be safe to compare the ratios of slopes. 



3.2 



Neves et al. ( 2009 ) Data 



Neves et al. (20091 published abundance data for 451 stars 



comparable in quality to ours. They measured the same el- 
ements we did (with the addition of Mg and Mn). We ap- 
plied the same culling criteria and analysis method as de- 
scribed above to their data, resulting in subsamples with 53 



the final corrected results in Figure 2 (see Gonzalez et al 



^ Note, we did not include the [Neves et al. abundance 
values for Mn in our analysis, as they tend to have larger uncer- 
tainties than the other elements they measured. 
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Figure 3. Average [X/Fe] values for the SWP (dots) and com- 
parison (open circles) solar analog stars from our samples. The 
average error bars show the standard error of the mean for the 
SWPs and comparison stars; they were calculated from the star- 
to-star scatter (upper left). 
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Figure 4. Same as Figure 2 but using the abundance data from 
the |Neves et al.||"2009| subsamples. 



SWPs and 225 comparison stars. We show the bias-corrected 
weighted-average abundance-Tc slope differences in Figure 
4. 

While there are fewer SWPs plotted in Figure 4 com- 
pared to Figure 2, similar negative average slopes are evi- 
dent in the figures. The average SWP [X/H]-Tc slope dif- 
ference value from these data is —2.3 x 10~^ ± 11 (s.d.) 1.5 
(s.e.m.) X 10~^ dex K"'^; 36 SWPs have negative slopes, and 
18 have positive slopes. 




Figure 5. Same as Figure 3 but using the abundance data from 
the [Neves et al.| | |2009[ l solar analogs. 



Three SWPs (HIP 15527 8033 7, 116906) and 14 com- 
parison stars from Neves et al. ( 2009 1 qualify as solar analogs 



according to our criteria. We plot the average relative abun- 
dances, as [X/Fe], against Tc for the SWP and comparison 
solar analog stars in Figure 5. 

The slopes of least-squares fits to the trends in Figure 
5 are (8.0 ±1.4) x lO"'^ and (15 ±4) x lO"'^ dex K"^ for the 
comparison and SWP solar analogs, respectively. This gives 
a ratio of 1.9 for SWPs relative to comparison stars, again 



consistent with the results of Melendez et al. ( 2009 1 



3.3 Combined Data 



There are 22 stars in common between Neves et al. ( 2009 1 



and the present work. We find good agreement for the Tefi 
and abundance-Tc slope values for these stars. The mean 
difference in Tog is — 16 ± 40 K. We compare the slopes in 
Figure 6. 

We adjusted our Teft and [X/H]-Tc slope values to be 
on the same scale as the Neves et al.| ( 2009 1 data based 



on the above comparisons. We then averaged the data for 
the stars in common and formed new combined SWP and 
comparison stars samples, which contain 100 and 277 stars, 
respectively. The bias-corrected combined [X/H]-Tc slope 
difference values are plotted in Figure 7. 

The average SWP abundance-Tc slope value from these 
data is -4.0 x lO"'^ ± 11 (s.d.) 1.1 (s.e.m.) x lO"'^ dex K'^ 
67 SWPs have negative slopes, and 34 have positive slopes. 
The number of stars with positive and negative slope val- 
ues are comparable below about 5800 K, but the negative 
slopes strongly dominate above 5800 K. The Sun's slope dif- 
ference is negative (—4.0 x 10~^ dex K~^), but stars with 
both positive and negative slopes are present near the solar 
temperature. 
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Figure 6. Comparison of the [X/H]-Tc slope values calculated 
from the abundances in the present work and from those of |Neves] 
|et al.| l |2009^ for stars in common. The solid line represents a one- 
to-one correlation; the dashed line is a least-squares fit. Typical 
error bars are shown in the lower right. 




Figure 8. Difference between the solar photospheric and mete- 
oritic elemental abundances is plotted against condensation tem- 
perature using the abundance data from Table 1 of |Asplund et alT] 
l |2009l l. The error bars include the errors from both the solar and 
meteoritic abundances. The solid line is a weighted least-squares 
fit. 
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Figure 7. Same as Figure 2 but using the combined samples. 



3.4 Solar System Abundances 



Gonzalez ( 1997 1 first noted a trend between the difference in 



the Solar System photospheric and meteoritic abundances 
and Tc . Gonzalez ( 2006 1 revisited this anomaly using the 



abundance data of Lodders ( 2003 1 and Asplund et al. ( 2005 1 



and found that both data sets display a significant trend. 
Since 2006, new photospheric and meteoritic abundance 
data have been published. We repeat the analysis of lGonzaT] 
[lei] ( [2006| ) below with more recent data to determine if the 
trend persists. 



We show in Figure 8 the difference between the solar 
photospheric and meteoritic abundances plotted against Tc 
using the data from Asplund et al. (20091. The figure in- 



cludes abundance differences with total uncertainties less 
than 0.10 dex. The slope from a weighted least-squares so- 
lution is (7.3 ± 2.5) X 10~^ dex K~^, and the Pearson's cor- 
relation coefficient is 0.43. The probability that the slope is 
actually zero is about 1%. 

We plot in Figure 9 the difference between the solar 
photospheric and meteoritic abundances using the data com- 



piled by Lodders et al. ( 2009 1. The best-fit slope to the abun- 



dance differences with total uncertainties less than 0.10 dex 
is (4.2 ± 1.8) X 10~^ dex K~^, and the Pearson's correlation 
coefficient is 0.37. The probability that the slope is actually 
zero is about 2%. 

There is reason to suspect that the data plotted in Fig- 
ure 9 underestimate the trend with Tc. The data compiled 
by [Lodders et al.| 12009] ) is a heterogeneous collection of re- 
cent abundance determinations from the literature. When 
multiple solar photospheric abundance determinations are 
available for a given element, [Lodders et al.[ ( |2009[ ) often 
select the values that are in closest agreement with the 
meteoritic value. For example, when discussing one of the 
published the Na photospheric abundances, [Lodders et al.| 
(2009]) write, "We do not adopt this value, as it is 25% 
lower than the meteoritic value as well as previously de- 
termined photospheric Na abundances." They make similar 
statements for K and Os. On the other hand, for P, they 
note that the new abundance value agrees well with the 



" The Co abundance we adopt in our analysis is from|Bergemann 
|et al.| | |2009| , which was published too late for |Lodders et al, 
1 2009 I to include it in their compilation. 



Trends with Condensation Temperature 5 




1000 1200 1400 1600 

T,(K) 



Figure 9. Same as Figure 6, but using the data from Table 4 of 
[Lodders et al.| l |2009| |. 



meteoritic value. This approach necessarily biases the abun- 
dances towards better agreement with the meteoritic values. 

If one takes a longer historical view of the study of So- 
lar System abundances, it is true that as the uncertainties 
have been reduced for many elements, the photospheric and 
meteoritic abundance values have come to agree much more 
closely. It is a mistake, however, to assume that the photo- 
spheric and meteoritic abundances will continue to agree as 
the abundance data continue to improve. 

Gonzalez (20061 obtained slope values of (8.2 ± 2.7) x 
10"^ dex K ' and (6.1 ± 2.5) x 10"^ dex K"^ from the 



Asplund et al. ( 2005 1 and Lodders ( 2003 1 data, respectively 



As a kind of control, we repeated the analyses shown 
in Figures 8 and 9 with atomic number, Z, replacing Tc. 
The resuhing slopes are (-3 ± 300) x 10" *^ dex Z"^ and 



(3 ±3) X 10""* dex Z~^ from the Asplund et al. (20091 and 



Lodders et al. 



( 2009 1 data, respectively. Neither of these 
trends is statistically significant. This gives us additional 
confidence that the trends in Figures 8 and 9 are real. 

Finally, we examine the solar abundances determined 
by [CafFau et al.| ( |2010| ). They measured only 12 elements, 
including Li, C, N, O. Excluding these four elements and 
another element with a large uncertainty (Os) , leaves us with 
the following seven elements: P, S, K, Fe, Eu, Hf and Th. A 
weighted least-squares fit gives a slope of (8.5 ± 6.6) x 10~^ 
dex K-\ 



4 DISCUSSION 

Our abundance data and the abundance data of INeves eti 



al. ( 2009 1 separately show that S WPs have more negative 



[X/H]-Tc slopes than stars without known planets; the com- 
bined data display similar patterns. These results imply that 
SWPs are relatively more depleted in high-Tc elements than 
low-Tc elements relative to comparison stars. 

[Melendez e\~aX\ ([2009| were the first to detect signifi- 
cant correlations between abundances and Tc among Sun- 
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Figure 10. Same as Figure 7, but only showing stars with [Fe/H] 
> 0.10. 



like stars. Restricting their study to solar analogs and twins, 
they found that the stars in their sample (with and without 
planets) have more positive slopes that the Sun. We con- 
firmed this finding with our examination of the solar analogs 
in our sample and the 



Neves et al. ( 2009 1 sample 



Melendez et al. (20091 also found that solar analog 



SWPs have more positive slopes than stars without planets. 
We confirmed this result in Figures 3 and 5, with an aver- 
age slope ratio between SWPs and comparison stars (~ 2.2) 
nearly identical to the ratio determined by [Melendez et aL] 
(20091. 



In a follow-up study, Ramirez et al. ( 2009 1 confirmed the 
findings of [Melendez et al.| ( |2009| ) with additional observa- 
tions of solar analogs and noted another pattern. They find 
that the [X/Fe]-Tc slope tends to be more negative among 
the more metal-rich stars. To search for this pattern in our 
data we have divided the data points from Figure 7 into two 
groups, [Fe/H] > 0.10 and [Fe/H] < 0.10, which we show in 
Figures 10 and 11, respectively. 

The differences between the trends in Figures 10 and 11 
are dramatic. In Figure 10 there are 48 SWPs with negative 
and 15 with positive slopes; restricting the count to values 
above 0.00007 and below -0.00007 dex K~\ there are 32 
SWPs with negative and 5 with positive slopes. In Figure 11 
there are 17 SWPs with negative slopes and 21 with positive 
slopes; restricting the count to values above 0.00007 and 
below -0.00007 dex K"\ there are 7 SWPs with negative and 
10 with positive slopes.. These results confirm fRamirez et al.| 
(2009) but for SWPs over a much more broad temperature 
range. 

From Figure 11, we can see that the Sun has a more 
negative slope than most SWPs near the solar temperature. 
Only one S WP in the figure (H IP 116906, one of the solar 
analogs from |Neves et al.] ( |2009[ )) has a negative slope (actu- 
ally very close to the solar value). Several SWPs hotter and 



" Note, the [X/Fe]-Tc slope values reported by [Melendez et al.| 
1 2009 I must be multiplied by —1 before comparing to our values. 
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Figure 11. Same as Figure 7, but only showing stars with [Fe/H] 
s£ 0.10. 



a few cooler than the Sun also have negative slopes. Only a 
larger sample over a similar range of temperatures will show 
if this is a statistical fluke or a significant anomaly for the 
Sun. 

Our new results for SWPs, combined with our recent 



results from Gonzalez (20081 and Gonzalez et al. (20101 



concerning Li abundances and vsini, likely share the same 



origin. As we suggested in Gonzalez (20081, low Li abun- 



dance and slow rotation among SWPs likely resulted from 
the presence of a massive protoplanetary disk around each 
star; Bouvier (20081 has also argued for such a relationship. 



While modern stellar evolution models that incorporate non- 
standard mixing are able to reproduce the low Li abundance 



of the Sun and SWPs (e.g., |Do Nascimento et al.| ( |2009[ )), 
it is still necessary to invoke another parameter to explain 
why single stars have different rotational histories. 



Ramirez et al. ( 2009 1 and Melendez et al. ( 2009 1 sug 



gest that refractory elements are relatively more depleted in 
the Sun because they were removed to form terrestrial plan- 
ets from the gas in the protoplanetary disk, gas which was 
later accreted onto it. Perhaps a similar process occurs for 
the SWPs. It is also important to note that hotter dwarfs 
have less massive envelopes. Therefore, hotter SWPs should 
be more sensitive to alteration of surface abundances by ac- 
cretion. There appears to be evidence for this in Figure 7. 

The correlation evident in Figures 8 and 9 are relevant 
to the interpretation of the results of Figure 11. It can help 
us to determine the source of the Sun's apparently anoma- 
lous composition. If refractory elements are relatively more 
depleted in the Sun because they were removed to form ter- 
restrial planets from gas that was later accreted onto it, then 
Figures 8 and 9 are telling us that some other process must 
also be at work. These figures show that the Sun's photo- 
sphere is enriched in refractory elements relative to the CI 
meteorites. If, instead the CI meteorites were depleted in re- 
fractory elements when they formed, the responsible process 
would have had to operate very early in the protoplanetary 
nebula, given the very primitive nature of these meteorites. 



In the case of the Solar System, however, the connection 
to terrestrial planets should be made using the meteorites 
formed in the inner Solar System (Alexander et al. 2001| 



and not the primitive CI meteorites that probably formed 
in the outer Solar System and at earlier times. Since the 
inner Solar System meteorites are enhanced in refractory el- 
ements by about a factor of two relative to the CI meteorites 
(Alexand er et alT]|2001 |, and since we have shown that the 
Sun is slightly (~ 0.05 dex) enhanced in refractory elements 
with respect to the CI meteorites, we conclude that the Sun 
is deficient in refractory elements (as already suggested by 
Melendez et al. (20091). The sequestration of refractory el- 
ements into the terrestrial planets was likely an important 
process. 

When discussing elemental abundance trends among 
stars, it is important to consider the effects of Galactic chem- 
ical evolution. Due to the design of our samples and analysis 
method, however, chemical evolution effects should not be 
a major infiuence on our results. This is because we have 
been careful to compare SWPs to otherwise very similar 
stars lacking Doppler detectable planets. Nevertheless, it is 
still possible that there is a mix of thick and thin disk stars 
in our samples, especially for the metal-poor stars. Thick 
disk stars have enhanced a-element abundances compared 
to thin disk stars for the same value of [Fe/H]; the possi- 
ble effects of this difference on planet formation has been 
discussed in Gonzalez (20091. How the different mix of ele- 
ments in a thick disk star afi'ects the [X/H]-Tc slope (e.g., 
for the data in Figure 11) is beyond the scope of the present 
work, but it is one that should be addressed as the sample 
size of metal-poor SWPs increases. 



5 CONCLUSIONS 

Using new abundance analyses of SWPs and stars without 
known planets, we have found that SWPs tend to have more 
negative [X/H]-Tc slopes than stars without planets. Our re- 
sults confirm [Ramirez et al.| ( [2C)09) , who focused their study 
on solar analogs. 

We also find that SWPs with [Fe/H] > 0.10 tend to 
have more negative [X/H]-Tc slopes than more metal-poor 
SWPs, showing that the process responsible for these trends 
is sensitive to metallicity. 

We revisited the question of the abundances in the so- 
lar photosphere relative to CI meteorites and confirmed pre- 
vious findings that a significant trend with Tc exists. The 
Sun is slightly (~ 0.05 dex) enhanced in refractory elements 
relative to the CI meteorites, but compared to the inner 
Solar System meteorites, the Sun is deficient by almost a 
factor of two. This implies that sequestration of the refrac- 
tory elements into terrestrial planets left their marks in the 
distribution of Solar System abundances. 

It appears that both the enrichment of refractory ele- 
ments in the solar photosphere via accretion and the seques- 
tration of refractory elements into the terrestrial planets left 
their marks in the distribution of Solar System elemental 
abundances. 

These results, combined with our recent findings that 
SWPs have lower Li abundances (confirmed by|Israelian et| 



al. ( 2009 1 for solar analogs) and rotate slower than compar- 
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ison stars, place new stringent constraints on planet forma- 
tion models ( Bouvier|2008 1. 
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